Fluid adsorption at a non-planar wall: 
roughness induced first-order wetting 
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We study the problem of fluid adsorption at a non-planar wall with a view to understanding the 
influence of surface roughness on the wetting transition. Starting from an appropriate Landau-type 
free energy functional we develop a linear response theory relating the free energy of the non-planar 
system to the correlation functions in its planar counterpart. Using this approach we are able to 
generalize the well known graphical construction method used to study the planar surface phase 
diagram and derive analytical expressions for the shift in the phase boundary for first and second- 
order wetting transitions. The results of the calculation are compared and contrasted with simple 
phenomenological and scaling arguments. Of particular interest is the influence of surface roughness 
on a second order wetting transition which is driven first order, even for small deviations from the 
plane. 
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While the statistical mechanical theory of fluid adsorp- 
tion at planar walls and in other idealized geometries 
(such as capillary-slits and cylindrical-pores) is much 
studied [QJ^] the microscopic theory of fluid adsorption 
at non-planar (rough) walls is far less developed. Nev- 
ertheless this problem is certainly of practical interest as 
well as posing something of a theoretical challenge due 
to the loss of translational invariance. For adsorption 
at a single wall perhaps the most important issue is 
whether the roughness influences the order and location 
(phase boundary) of any wetting transition ||. Here 
we address this question using a generalized Landau 
density-functional model which may be viewed as the 
simplest available microscopic approach [Q. Recall that 
in application to theories of the planar liquid-vapour 
interface as well as wetting transitions and finite-size 
effects the same approach has played a pivotal role jj| 
leading to more sophisticated methods. It is therefore a 
natural starting point for the systematic investigation of 
non-planar fluid interfaces. A major part of our work is 
to show that within perturbation theory (valid for small 
deviations from the plane) it is possible to derive an 
analytical expression for the change in free energy due 
to roughness that may be studied using a generalization 
of a graphical construction method familiar from the 
planar problem In this way we avoid having to 

use effective Hamiltonian methods [Q. Our analysis 
borrows results and methods recently developed for the 
calculation of pair correlation functions at planar wall- 
fluid interfaces ||] and is we believe, of pedagogical 
interest beyind its application reported here. Using this 
approach we are able to derive analytical expressions for 
a shift of the phase boundary at first and second-order 
wetting transitions in three dimensional systems. While 
our analysis is mean field-like we are confident that the 
topology of the surface phase diagram is uneffected by 
the inclusion of fluctuation effects. 

To begin we make some preliminary remarks con- 
cerning the possibility of a roughness induced wetting 



temperature shift. It was observed by Wenzel |)| some 
sixty years ago that the contact angle 9^ of a drop of 
planar liquid (phase (3 say) on a rough substrate (wall-o 
phase interface) appeared to satisfy the empirical relation 



cos 9-e — t cos 0^ 



(1) 



where r = 4 s is the ratio of non-planar to planar surface 
area and 9 V is the contact angle for the planar system 
(r = 1). Of course this is not an exact result but has 
nevertheless been proved (rigorously) valid for an Ising 
model with a non-planar boundary at sufficiently low 
temperatures [ [To| |. Unfortunately the low temperature 
restriction precludes a study of the influence of roughness 
on any wetting transition. However it is clear that 
if Wenzel's law was valid V6>h > then the wetting 
temperature is necessarily reduced. This is conveniently 
expressed as follows : the non-planar wall-a interface 
is completely wet (i. e. 9s = 0) by a fluid /3 phase for 
roughness parameter r > ryy where 



r w 



sec 9„ 



(2) 



and recall 9 V is the contact angle (> 0) of the (3 
droplet at the planar wall-a interface (at two phase bulk 
coexistence) . If we assume that the wall has a corrugated 
shape described by the graph zjy( r j|) = V^asinqx where 
r 1 1 = (x, y) is the displacement vector parallel to the z = 
plane (corresponding to the mean position of the wall) 
then we can rewrite the equation for the phase boundary 
as 



qa 



(3) 



assuming 9^ is small. Here we have included a factor \J~2 
so that a measures the root mean square width of the 
wall. We shall refer to equations (^,||) as Wenzel's result 
although of course Wenzel was unaware of the possibility 
of a wetting phase transition. Thus it appears to be 
possible to induce wetting by increasing the roughness 
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of the substrate although this becomes increasingly more 
difficult for planar contact angles close to f . However 
one should be suspicious of this prediction given that 
this approach makes no mention of the order of the 
wetting transition occuring in the planar system (at 
temperature T v say). In particular, experience with the 
well developed finite-size scaling theory of bulk critical 
phenomena suggests that a more reliable expression for 
rw would contain information about the (transverse) 
correlation length, which recall diverges at a second order 
wetting transition. Similarly at this simple level we are 
not able to offer any prediction whether the order of the 
wetting transition in the non-planar system is different to 
that occuring for planar geometry. In fact, as we shall see 
the Wenzel result is inaccurate as regards the influence 
of roughness on second-order wetting transitions which 
turns out to be much more interesting than ([|) suggests. 
However, we are able to provide a microscopic derivation 
of (0) and(||) for first order wetting. 

Before we discuss the microscopic Landau theory we 
note that it is straight forward to develop a scaling theory 
for the roughness induced phase shift which improves 
on (^). Let us suppose that in the planar system the 
wall-a interface undergoes a continuous (second-order) 
wetting transition as t = T ^T T — > 0. By definition the 

2-a g 

contact angle 9^ vanishes as 9 n ~ t = where as is the 
specific heat exponent characterising the singular part 
of the excess free energy / snlg ~ t 2 ~ as . Assuming that 
the wall has the corrugated shape mentioned above it is 
natural to expect that the length scales a and q enter 
the singular part of the non-planar free energy through 
scaled variables at^ s and qt~ v ^ where (3s and v\\ are 
the adsorption and transverse correlation length critical 
exponents respectively j| . In this way we deduce that the 
phase boundary in the non-planar geometry (implicitly 
incorporating the wetting temperature shift) satisfies the 
scaling relation 

T (t,...)~ q *i L S{<un i ) (4) 

where we have used the hyperscaling relations 2 — as — 
(d— l)f|| and (3s — ^rV\\ valid below (and at) the upper 
critical dimension d* (which recall is equal to three for 
systems with short range forces jjj). Here S(w) is an 
unknown scaling function which must satisfy S(0) = 
in order to reproduce the trivial planar phase boundary 
Ok = 0. Interestingly Wenzel's result (||) does meet 
with the scaling requirement (§) although as we shall 
see the correct mean-field scaling function S(w) is more 
complicated and shows non-analytic behaviour. With 
these preliminaries in mind we now turn to the analysis 
of a microscopic density-functional model. 

Consider a one component system with order param- 
eter m(r) which shows bulk two-phase coexistence at 
sub-critical temperatures (T < Tc) in zero bulk field 
h = between phases with order parameters m a (T) (> 
0) and mp(T) (< 0) respectively. If the system is 
bounded by a fixed wall whose position is specified 



by a height variable zvK r ||) we suppose that the free 
energy functional F[m(r)] accounting for bulk and wall 
interactions is 

F[m(r)] = J dr j^(Vm) 2 + <j>(m) 

+S(z - z w (r n ))[l + i(Vz w ) 2 ]0i(m)| (5) 

which naturally generalizes the standard free energy 
model of the planar semi-infinite system @j. As usual 
the bulk free energy term <j)(m) has a double well form 
for T < Tc but will not be specified further. The surface 
interaction term 4>i {m) p is taken to have the standard 

expression <f>i(m) = him where c is the surface 

enhancement and hi is the surface field. Note that the 
surface term contains an extra factor 1 + i(V zw) 2 which 
accounts for the increase in surface area due to deviations 
in the position of the wall from the plane 0. We 
will always assume that these deviations are small and 
also of long wavelength compared to some appropriate 
microscopic scale (see later). 

It is natural to anticipate that for small deviations from 
the plane the minimum of _F[m(r)] may be written as a 
perturbation about the planar value 

i*E = <t>{m b )V + oA^ 

+ 2{2 l )d -i J dc L <7 2 A.(g)|Mq)| 2 + ' ' ' (6) 

where the ellipses denote higher order products of the 
Fourier amplitudes Zw{q) which can be safely ignored 
for small deviations. The first two terms represent the 
free energy of the planar system so that 4>(rrib) and a 
correspond to the bulk free energy density and surface 
tension respectively. The latter quantity is given by the 
well known expression [0J|^,§| 

/■"if, 

<r=0i(mi)± / dmQ (m) (7) 

J mi 

where m& and mi are bulk and surface magnetizations 
respectively, and 

Qo(m) = y/2[<t>(m) - <Mm 6 )] (8) 

The sign in (j7|) is chosen so that the contribution from 
the integral is positive. This expression is amenable to 
a well known graphical interpretation which is extremely 
useful in determining the (planar) surface phase diagram 

1- 

The quantity A^g) to be determined represents the 
free energy correction due to roughness and has the 
same dimensions as the surface tension a. Note that 
if A„.(<7) were wave vector independent then the free 
energy increment is simply proportional to the increase 
in surface area — A n . 

The starting point in our analysis is an exact linear 
response relation for A w (q) in terms of the (Fourier trans- 
formed) planar pair correlation function G(i"i,r 2 ) =< 
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m(ri)m(r2) > — < m(ri) >< mfa) > where both 
particles are exactly at the wall: 



G(0,0;q) = / dry G(r u r 2 



(9) 



with ri = (0,0) and r 2 
conveniently written 

q 2 A n (q) = g 2 </>i(mi) + m 



(0,r||) 



The result is 



G(0,0;q) 
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(10) 



where m'i is the derivative of the (planar) wall magnetiza- 
tion with respect to z and xu is the surface susceptibility 
Xn = While a derivation of ( |Io| ) is not appropriate 

here some remarks should clarify its origin. The 
first term is the 'direct' potential contribution due to 
the increase in surface area. Note also that the term 
in parentheses vanishes when q = (by virtue of a 
sum rule) reflecting the arbitrariness in the location of 
the z = plane. The all important dependence on 
G(0, 0; q) follows from remarks made earlier by one of us 
Jl2| concerning the exact relationship between G(z : z; q) 
and the free energy cost of a fluctuation in the location 
of a surface of (appropriately) fixed magnetization m x 
whose average position is z. Fortunately a good deal 
of information about G(0, 0; q) at wetting transitions 
is now known so that ([LG) constitutes a rather useful 
relation. In fact it is possible to continue further 
with the analysis for arbitrary <f>(m) (and 4>i(m)) and 
derive an elegant expression for A^^q) comparable 
with the surface tension formula (Q). To proceed we 
substitute into (|io| ) the known, exact expressions |ljj 
for the moments G2 n (0,0) appearing in the expansion 
G(0,0;q) = J2n=o 9 2 ™ G 2n(0, 0). While these moments 
are themselves rather cumbersome the corresponding 
moments A 2 „ in the expansion A„-(g) = Y^=o 1 2n 
turn out to be much more compact. In this way we 
we arrive at the following equation for A^g) valid for 
arbitrary 0(m) and </>i(m) 



Aiv(q) = <M m i) ± / dmQ(m;mi,q) 



(11) 



where the sign is the same as that in (|7j). Here the new 
function Q(m;mi,q) satisfies the integral equation 

Q(m;mi ,q) 

Qo(m) 

1 - q 2 Q dm'Qo V) dm"Q(m"; m ls q) (12) 

where we have assumed that > iiij so that the sign of 
the double integral is positive. The function Q(m; mi, q) 
satisfies the boundary conditions Q(mi,;mi,q) = 
and Q{mi]mi,q) — Qo(mi) which are useful when 
interpretating the equation graphically. An important 
result immediately follows from ([?]) , ( J]~l[ ) and ([l2]) namely 



Q„(m) 




X(m) 



FIG. 1. Illustration of the generalized graphical con- 
struction method for first-order wetting. The dashed lines 
represent the non-planar Q functions while the shaded areas 
show the contributions to the free energy from the second 
term in (11) for wet (W) and partially wet (P) profiles. 



This identity reflects the invariance of the free energy 
of a planar wall-fluid interface with respect to rotations 
and is related to the asymptotic coherence of surface 
correlations discussed at length in We have also 

found it profitable |llj to consider the differential 
version of ( [l2] ) conveniently expressed in terms of the 

dimensionless scaling factor y(m;mi,q) = — q^'^' 9 ^ 
which satisfies the boundary conditions y(m\\ mi, q) = 1 
and y (m&; mi , q) = 



Ql{m)y" + 3Q (m)Q' (m)y' -q 2 y = 



(14) 



where prime denotes differentiation with respect to m. 
For the standard l (/> 4 ' theory (in zero bulk field h = 0) 
it transpires that y has a scaling structure y(m; mi, q) = 
y(9«b)(JTL ; ™i) where £& = k^ 1 is the bulk correlation 

length. The function Y^(t;ti) satisfies 



(1 - t 2 )Y^ - 6tY (a *> ~ 



y(«) = o 



(15) 



A,r(0) = a 



(13) 



which is an adaptation of Latta's generalized Mathieu 
equation Jh|. Here the overdot implies differentiation 
with respect to t. Equations ( |l0| - |l5| ) constitute the main 
results of our perturbation theory which we now apply 
to the problem of wetting at a non-planar wall. 

An important result follows directly from ([h]). As- 
suming that the wetting transition occuring in the 
planar geometry is second-order it follows that the 
correlation function G(0, 0; q) and hence the non-planar 
free energy cannot exhibit singular behaviour away from 
T n . Consequently the wetting transition is shifted only if 
the transition is roughness induced first-order. Similarly 
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a first-order wetting transition in the planar geometry 
cannot be roughness induced second-order since this 
would require an unphysical singularity in G(0,0;g). 
With these preliminary remarks in mind we consider the 
case of the influence of roughness on strongly first-order 
and second-order wetting transitions seperately. 

a. Strongly first-order In figure 1 we illustrate 
a useful graphical interpretation of ( jL2| ) for A 7r (g). 
Intersections of the straight line X(m) — cm — hi with 
Qo(m) determine the surface magnetizations ro^mf of 
the wet and partially wet planar magnetization profiles 
respectively. Thus the non-wet profile starts at mf and 
increases to m a while the wet profile starts at m^ and 
increases to mp (with a (3a interface at infinity). 

The function y(rn;mi,q) can be found accordingly 
for these profiles. It is very nearly unity (provided 
q£b <C I) for the non-wet profile and also for the wet 
profile between m^ and rap while it is zero between 
To/5 and m a . The corresponding Q functions are also 
illustrated in figure 1 (in an obvious notation) together 
with their respective integral contribution to A n (q) 
shown as shaded areas. In this way we conclude that the 
correction term A p for the partially wet profile is ~ 
a wa while for the completely wet profile it is ~ cr w p. 
Note that the surface tensions of these profiles are given 
by a p — a wa and a w = a w p + a a p respectively. Then 
a simple free energy balancing argument shows that the 
phase boundary is shifted exactly according to Wenzel's 
result (|) 

b. Second-order wetting For this case the analysis is 
more involved due to the presence of long wavelength 
fluctuations in the adsorbed fluid film. Consequently we 
only quote our main results which, for simplicity are 
restricted to the case where the wall is taken to be at 
position z\y = ^/2asinqx (and recall q£,b <C I). For 
fixed 9-x ~ T-x — T the planar second-order transition 
is roughness induced first-order for roughness parameter 
r satisfying 



rw ~ sec + (sec 9^ — l)q 2 



(16) 



where £y is the planar value of the transverse correlation 
length at this temperature. Clearly this is a fluctuation 
modified version of Wenzel's result (||). 

Note the factor > is a measure of the deviation 

C+K 

from the planar tricritical point. For small values of 9^ 
i. e. T close to T n the result reduces to 



r w ~ sec 6^ + R (q£ b ) 



where 



f? — lim 



/sing£j| 



(17) 



(18) 



is the ratio of two hyperscaling amplitudes and takes the 
universal value 72 = i in mean field theory. Thus the 




FIG. 2. Schematic wetting phase diagram for fluid adsorp- 
tion in a system with a non-planar wall. The lines C\ and 
C2 are loci of first and second-order wetting phase transitions 
respectively which meet at a tricritical point corresponding 
to a — a* < The vertical axis is a linear measure of the 
temperature scale T n — T for small 8 n . 



shifted phase boundary (for 9 V > 0) is given by the non- 
analytic function 



for a < a* 

q(a 2 — a* 2 )^ for a > a* 



where 



(19) 



(20) 



is the tricritical value of a. The analytical expression 
is confirmed by numerical minimization of the mean 
field free energy functional Thus for fixed a < 

a* the wetting phase transition remains second-order 
and occurs at the (unshifted) temperature T„. While 
for a > a* the transition is roughness induced first- 
order and occurs at a lower temperature. The phase 
diagram plotted in terms of 9^ and a is sketched in 
figure 2. Importantly, the expression for the shifted phase 
boundary is precisely of the scaling form (|j) provided 
we identify d — 3 corresponding to the upper critical 
dimension for wetting (for models with short ranged 
forces) as we might anticipate for the present mean field 
calculation. Interestingly as c is reduced to k (i. e. as 
we approach the planar tricritical point) we recover the 
simple prediction of Wenzel (||). It is also noteworthy 
that the tricritical value of the width parameter a* given 
by ( p0| ) is independent of q (for q£t *C 1). The fact that 
a second order wetting transition is roughness induced 
first-order for even minor deviations from the plane is 
the central conclusion of our study. Moreover experience 
with (improved) effective Hamiltonian models of wetting 
in three dimensional systems is strongly suggestive that 
this result is true beyond mean field theory since there 
does not appear to be a mechanism (involving a position 
dependent stiffness or stiffness matrix [^5|) by which 
a first-order wetting transition is fluctuation induced 
second-order. Thus we are confident that the topology of 
the surface phase diagram (see figure 2) would be borne 
out in Ising model simulation studies. 



4 



In summary, we have developed a mean field linear 
response theory for fluid adsorption at a non-planar wall 
and derived an exact analytical expression for the pertur- 
bation of the planar free energy due to surface roughness. 
Using this method we have been able to vindicate a 
phenomcnological result (which may be attributed to 
Wenzel) regarding the influence of roughness on first- 
order wetting transitions. For second-order wetting 
transitions however, Wenzel's result is inappropriate (due 
to fluctuation effects) and calculation shows that the 
transition is driven first-order for deviations from the 
plane. 
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